Intrinsic Structural Disorder and the Magnetic Ground State in Bulk EuTiOs 
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The magnetic properties of single-crystal EuTiC>3 are suggestive of nanoscale disorder below 
its cubic-tetragonal phase transition. We demonstrate that electric field cooling acts to restore 
monocrystallinity, thus confirming that emergent structural disorder is an intrinsic low-temperature 
property of this material. Using torque magnetometry, we deduce that tetragonal EuTiC>3 enters 
an easy-axis antiferromagnetic phase at 5.6 K, with a first-order transition to an easy-plane ground 
state below 3 K. Our data is reproduced by a 3D anisotropic Heisenberg spin model. 

PACS numbers: 75.85.-t-t, 75.30.Gw, 64.60.Cn, 75.10.Dg 



Perovskite titanates AT1O3 are among the brightest 
stars in the rapidly-developing field of oxide electron- 
ics. Their dielectric and transport properties are eas- 
ily modulated by epitaxial strain, field-effect or chemi- 
cal doping pH!] ; hence facilitating their integration with 
conventional charge-based electronics. However, the past 
decade has also seen an intensive search for a multiferroic 
perovskite with strong magnetoelectric coupling, which 
would enjoy numerous spintronic applications. The spot- 
light soon fell on EuTiC>3, whose large magnetic mo- 
ment (S — 7/2 per unit cell) and quantum paraelec- 
tricity suggest that it is a magnetic analogue of SrTiC>3, 
the "workhorse" oxide for electronic devices. Although 
EuTiC>3 has been known as a G-type antiferromagnet 
below Tat ~ 5.5 K since the 1960s [S], its magneto- 
electric properties were only revealed in 2001 by a 7% 
drop in the dielectric constant at Tjy JB]- Early ef- 
forts to model EuTiC>3 shared one critical feature: an 
assumption of cubic crystal symmetry throughout the 
phase diagram [?HS]. The recent discovery of a cubic- 
tetragonal transition at 283 K in powdered samples [10] 
obliges us to revise our views, since tetragonal symme- 
try pushes EuTiC>3 closer to a ferromagnetic/ferroelectric 
(FM/FE) phase boundary [7] and permits the breakage 
of spatial inversion symmetry. Theoretical [11 and X-ray 
diffraction (XRD) analysis [12] suggests that tetragonal 
EuTiC>3 may be disordered, thus opening the door to 
local symmetry-breaking and electric dipole formation 
via the Dzyaloshinskii-Moriya (DM) interaction. More 
generally, nanoscale heterogeneity limits applications for 
many complex oxides [13] . since functional electronic 
materials must display phase purity over lengthscales 
greater or equal to their intended device dimensions. It 
is therefore an urgent priority to identify methods for 
controlling or suppressing such disorder. In this Let- 



ter, we employ the magnetic anisotropy of EuTiC>3 as 
a structural probe, hence showing that monocrystalline 
cubic EuTiOs undergoes a transition to an intrinsically 
nano-disordered tetragonal phase with an easy-plane AF 
ground state. However, the disorder may be reduced by 
cooling under an electric field. This mechanism of homo- 
geneity control may prove crucial for phase management 
in future devices made from elastically hard oxides. 

Our experiments were performed on EuTiC>3 single 
crystals, grown in a floating-zone furnace [14j . High- 
symmetry axes were determined by XRD in the cu- 
bic phase at 300 K and the crystals cut into cuboids 
with faces parallel to (100), (010) and (001). The cut- 
ting/alignment error is ~ ± 5°, a figure supported by 
our torque data. We stress that at 300 K, our samples 
are perfectly monocrystalline with no evidence for twin- 
ning. Geometric demagnetization corrections [15j have 
been applied to all our data, except the specific heat and 
torque. Our results have been verified using four separate 
crystals, all of which display similar behaviour. 

The structural transition in EuTiOs is visible as a kink 
in the heat capacity C(T) at T c = 283 K (Fig. [IJi). As 
we sweep the temperature below T c , repeated thermal 
relaxations yield two divergent traces in C(T), the lower 
of which is generated by the first relaxation at a given 
temperature and vanishes below ~ 240 K. This trend is 
reproducible for all our crystals, but absent from refer- 
ence materials tested on the same calorimeter and hence 
not an experimental artifact. Such behaviour is sugges- 
tive of phase metastability in this temperature range. In- 
terestingly, XRD indicates a lower T c of ~ 235 K [12] . 
coinciding with the disappearance of our lower trace. 
We extract the temperature-integrated phonon density 
of states (PDOS) from C(T) and note the broad spectral 
weight distribution below 20 meV: this corresponds to the 
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FIG. 1: (a) Zero-field heat capacity C(T) in EuTi0 3 . Top 
left: zoom on the structural transition at 283 K, showing 
two divergent traces generated by consecutive thermal relax- 
ations. Lower right: PDOS extracted from the specific heat 
(see ref. [16] for details), (b) Low temperature magnon heat 
capacity C mag , illustrating the field-suppression of AF order 
and spin gap formation. The absence of the spin-flop and spin 
realignment transitions from C(T,H) is unsurprising, since 
the entropy change is extremely small and short-pulse relax- 
ation calorimetry is notoriously insensitive to first-order tran- 
sitions [17| . (c) dc magnetic susceptibility x(T) for a range of 
applied fields, revealing the transition to afe-plane order below 
~ 3 K, a "shoulder" due to the spin-flop and a peak at Tjv- A 
faint anomaly is also visible in low fields around 4.5 K (black 
arrow), (d) dc magnetization M(H) along the crystal axes. 
Small anisotropics in m aa t arise from an anisotropic g-tensor 
in a tetragonal system; the fact that our observed anisotropy 
is quasi-3D rather than 2D is indicative of disorder, (e) Spin- 
flop transitions in x(H) for fields parallel to each axis. In the 
absence of a spin-flop, ^ should remain constant down to 
H = (as illustrated by a dashed line for H // [010]). 



TOl phonon softening from 120 cm -1 to 80 cm -1 |18j . 
Subtracting a /?iT 3 + P2T 5 phonon background from 
C(T) at low temperature, we obtain the magnon heat 
capacity C mag (Fig. [T]b) . A large jump at T N — 5.6 K 
indicates a second-order AF transition, similar to exist- 
ing powder [TU] and nanoparticle [115] data. Integrat- 
ing Cmag/T from zero to 5.6 K yields a total entropy of 
69.4 ± 1 mJgat _1 K _1 which agrees with the theoretical 
value Nk B la(2S + 1) = 69.8 mJgat^K- 1 and discour- 
ages the possibility of any further structural transition 
at T < 5.6 K. Magnetic fields suppress Tjy, inducing 
a T = transition (at H = H sat ~ IT) to saturated 
paramagnetism, with an energy gap for spin excitations 
growing linearly in H — H sat . This confirms the phase 
purity of our samples, since any oxygen deficiency should 
drive a zero-field crossover from AF to FM order [20] . 

Further complexity emerges in the dc magnetic sus- 
ceptibility x — M/H (Fig. [lp) , where three energy scales 
may be identified. The peak in x(^) a t is robust at 
low fields; however a "shoulder" forms below T/v, broad- 
ening and moving to lower temperature as H increases. 
A third anomaly is visible below 3 K, vanishing above 
H ~ 0.1 T. The magnetization saturates at a similar 
field H sat ~ 1 T along all three axes (Fig. [TJi) , with a 
moment m sa t = 7.08 ± 0.2 /is/unit cell. At low fields, 
small anomalies in M{H) are visible and x{H) falls as 
H tends to zero (Fig. [TJj). The temperature and field- 
dependence of these non-linearities in M(H) correlate 
with the "shoulder" in x(T), leading us to identify them 
as spin-flop transitions. However, a single-crystal AF 
with tetragonal symmetry should only exhibit a spin-flop 
when H is parallel to the "easy" axis (and no spin-flop 
should occur in an easy-plane system). In contrast, our 
crystals consistently display spin-flops at similar fields 
(H sf = 0.22 ± 0.04 T at 2 K) along all three axes . 

To resolve this conflict, let us suppose that EuTiC>3 is 
indeed disordered below T c , as suggested by our C(T) 
and M(H) data as well as other recent work [TTJ[TJ]. I n 
this case, tetragonal EuTiC>3 may be considered as effec- 
tively polycrystalline. Our experiments do not permit us 
to estimate the lengthscale of this disorder; however, we 
note that spontaneous twinning is observed in elastically 
soft SrTiC>3 with domain size below 1 /zm |21j . The pres- 
ence of a spin-flop in a 200 nm EuTiC>3 film also suggests 
that true nanoscale disorder exists in this material |22j . 
Nevertheless, it should be possible to reduce this disorder 
by cooling EuTiC>3 with an electric field E applied par- 
allel to a cubic crystal axis: this prevents softening for 
phonon wavevectors k//E [23] and hence ensures that 
the tetragonal c-axis is perpendicular to E. It is crucial 
to -E-field-cool (EFC) from above 50 K, since below this 
temperature the phonon softening terminates, quantum 
fluctuations dominate the low-energy phonon spectrum 
and the lattice structure is presumably locked. Fig. [2] 
shows the effect of EFC on x(H) for a EuTiC-3 crystal 
cooled from 80 K: the spin-flop is clearly suppressed, im- 
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FIG. 2: Electric field-induced spin-flop suppression in an 
EuTiO-3 crystal cooled from T = 80 K under E = 230 V cm -1 
(E.LH). To apply the field, gold contacts were sputtered on 
opposing crystal faces and wires attached with silver epoxy. 

plying that an increased proportion of our crystal now 
exhibits c//H. This confirms the presence of structural 
disorder and its influence on the magnetic ground state. 

Since europium compounds exhibit significant neu- 
tron absorption and are hence poorly-suited to neutron 
diffraction, torque magnetometry is our best alternative 
to elucidate the spin structure [MJ |5S]. An easy-axis 
two-sublattice AF (Fig. [3^,) will generate an out-of-plane 
torque signal t{4>) oc sin20 when a field H rotates in a 
plane containing the easy axis; however r = if H lies 
in the plane perpendicular to the easy axis. Conversely, 
an easy-plane AF (Fig. [3]b) will display t(0) cx sin4# 
when H is rotated in the easy plane and t(4>) cx sin2</> 
when H is rotated through any plane containing the hard 
axis {4> and 9 are the polar and azimuthal angles respec- 
tively). Anisotropic paramagnetism above 7jv will also 
generate a t(</>) cx sin20 signal. 

Our magnetization data imply that EuTiC>3 crystals 
contain AF domains with the easy-axis pointing along 
all three crystallographic axes. Therefore, the torque 
should behave as r(a) = Asin2a + i?sin4a, where 
a is the field angle in a plane defined by any pair of crys- 
tal axes. (The demagnetization factor also contributes to 
t(o) with 180° periodicity, but it is temperature invari- 
ant and does not affect our conclusions.) A series of r(a) 
curves for 2 K > T > 15 K is shown in Fig. |3j;, with 
H = 0.05 T rotating in the (001) plane: the transition 
from sin2a behaviour above Tjy to sin4a at low temper- 
ature is clearly visible. Fitting our r(a) curves to obtain 
the sin2a and sin4cv weight parameters A and B, we note 
that B only appears below T/v, rising rapidly below 3 K 
with a corresponding drop in A. The only mechanism 
which can generate this behaviour is a change from easy- 
axis to easy-plane symmetry: we deduce that a first-order 
phase transition from c-axis to a&-plane AF ordering oc- 
curs close to 2.75 K. Upon increasing H (Fig. [3J3) , r(a) 
changes shape from sinusoidal to "sawtooth" , indicating 
a spin- flop [25!. At H — 0.25 T, the approximate sin2cv 
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FIG. 3: (a,b) Torque experimental schematic, with easy axes 
shown in red, AF sublattice moments in green and field H in 
blue, (c) Temperature dependence of t(ol) / / [001] for H = 
0.05 T. Rotating H in the (100) and (010) planes gives similar 
results, (d) Temperature dependence of the separated sin2a 
and sin4« contributions to the total torque shown in (c). (e) 
Field-evolution of the effective susceptibility x( a ) = T(a)/H 2 
revealing the spin-flop transition. 



dependence is characteristic of a spin-flopped AF. 

We summarise our low-temperature data in Fig. [4^,. 
For zero-EFC EuTiOa , little or no anisotropy is observed 
for different magnetic field orientations, as expected for 
an effective polycrystal. This phase diagram indicates 
that single-ion anisotropy terms induced by the tetrago- 
nal crystal field play an important role at low tempera- 
ture. To model the effect of these terms, we consider a 
classical Heisenberg model with nearest-neighbour (Ji) 



4 





1.2- 




1.0- 






~o 


0.8- 


03 
U- 


0.6- 


O 




nel 


0.4- 










5 


0.2- 




0.0- 



(a) 



-r Anisotropic 
* J t Paramagnet 

Spin-flopped r<H ^ 

Antiferromagnet |©| 
H 



ab-AF 



1^4^* "aF*** Hfljj 



(b)H//[100] 



Anisotropic 
Paramagnet 



* ♦ 4 

Spin-flopped * $ 



b-AF Antiferromagnet I 



+ + + 
C-AF 

-"i 1 1— 



(c)H//[001] 



Spin-flopped 
Antiferromagnet II 



Anisotropic 
Paramagnet 



C-AF 

t 



2 3 
Temperature (K) 



FIG. 4: (a) Experimentally-determined magnetic phase diagram of EuTiOs. (b,c) Two-sublattice mean-field phase diagrams 
for H II [100] and [001] (see text for parameters). afc-AF refers to zero-field AF order along the [110] or [110] axes, while c-AF 
describes AF order along [001]. For the spin-flopped AF I phase, the spins order in the ab plane with the staggered component 
of the magnetization perpendicular to the applied field. In contrast, spin-flopped AF II corresponds to AF ordering along 
[110] or [110], with a uniform magnetization component along the c-axis. The intrinsic structural disorder implies that our 
experimentally-determined diagram (a) should be a superposition of (b) and (c). Experimental evidence for the small pocket 
of c-AF in (c) is provided by the faint anomaly visible at low fields in x(T) (indicated by an arrow in Fig. [TJi.) 



and next-nearest-neighbour (J 2 ) exchange interactions 

i~l = j\ S r • S r ' -\- J2 ^ S r • s r > 

(r,r') {{ r j r ')) 

+ J2v(s r )-J2B-S r , (1) 

r r 

together with a single-ion anisotropy potential 

V(s) = -C 2 (S 2 + S 2 y )-C i (Sj + S*), (2) 

where S x = S ■ (1, 1, 0)/-y/2, and S y = S ■ (-1, 1, 0)/\/2. 
S r represents classical spin at r with IS^-I = 7/2. The 
exchange parameters were determined by fitting the mea- 
sured T — saturation field H sat — 12SJi/g/j,B ~ 1 T 
and Curie- Weiss temperature Tew = — 2S 2 ( J\ + 2J 2 ) ~ 
4 K [S] (g ~ 2 is the gyromagnetic factor), resulting in 
J1S 2 = 0.31 K and J 2 S 2 = -1.22 K. This yields a mean 
field T N = 2S 2 {Ji - 2J 2 ) ~ 5.5 K, in excellent agree- 
ment with the measured value. The single-ion anisotropy 
parameters C 2 S 2 = -0.4 K and C4S" 4 = 0.455 K were 
chosen to reproduce our experimental phase topology, es- 
pecially the transition at T a b ~ 3 K. The phase diagrams 
obtained from a two-sublattice mean-field theory with 
magnetic field applied along the [100] and [001] directions 
are shown in Figs. [4]b and c. Although our calculated c- 
AF region persists to slightly higher fields than our data, 
the overall agreement is reasonably good considering the 
limitations of a mean field approach. 

To the best of our knowledge, EuTiOs is the first ex- 
ample of a material undergoing a transition from a crys- 
talline to disordered state at low temperature, despite 
lacking any intrinsic structural, magnetic or electronic 
frustration. It is natural to question the effect of this 
disorder on the dielectric properties of EuTiOs: although 
neither the structural transition nor the AF order alone 
provide the spatial inversion symmetry breakage neces- 
sary to stabilise ferroelectricity, it is plausible that local 
disorder could cant the Eu 2+ moments via the DM inter- 
action. This would drive a local symmetry-breakage and 



hence enable the formation of nanoscale electric dipolcs. 
Structural disorder implies that the global dipole mo- 
ment would average to zero, rendering FE order difficult 
to detect. However, the possibility of a static polariza- 
tion extending over useful lengthscales in less-disordered 
EuTiC>3 merits further attention. 
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